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Abstract. We present herein the most recent BABAR results on direct CP asymmetry 
measurements in B — > X s y, on partial branching fraction and CP asymmetry measurements in 
B —» X s l + l~ , on a search for B —> n/r/£ + £~ decays, on a search for lepton number violation in 
B + —>• X~£ + £' + modes and a study of B° —> ujuj and B° —» uj(j> decays. 


1. Introduction 

The decays B —> X s ^7 and B —> X s ^£ + £~ (with £ + £~ = e + e _ , R + p~) are flavor-changing 
neutral-current (FCNC) processes that are forbidden in the Standard Model (SM) at tree 
level. They occur in higher-order processes and are described by an effective Hamiltonian 
that factorizes short-distance contributions in terms of scale-dependent Wilson coefficients 
CM 0 front long-distance contributions expressed by local four-fermion operators Oi that 
define hadronic matrix elements, 


H eS = 



Ci{p)Oi. 


( 1 ) 


While Wilson coefficients are calculable perturbatively, the calculation of the hadronic matrix 
elements requires non-perturbative methods such as the heavy quark expansion (2jjl]. 



Figure 1. Lowest-order diagrams for B —> X s ^-y (left) and B —>• X s ^£ + £ (middle, right). 


Figure [l] shows the lowest order diagrams for these FCNC decays. In B —> -W,d7) the 
electromagnetic penguin loop dominates. The short-distance part is expressed by the effective 
Wilson coefficient Cy®. Through operator mixing at higher orders, the chromomagnetic penguin 
enters whose short distance part is parameterized by Cf®. In B —> X Si rf£ + £ _ modes, the 
Z penguin and the WW box diagram contribute in addition whose short-distance parts are 
parametrized in terms of Cg® and Cf®, the vector and axial-vector current contributions of 
these diagrams. Physics beyond the SM introduces new loops and box diagrams with new 
particles (e.g. a charged Higgs boson or supersymmetric particles) as shown in Fig. [ 2 ] (left, 
middle). Such contributions modify the Wilson coefficients and may introduce new diagrams 
with scalar and pseudoscalar current interactions and in turn new Wilson coefficients, Cs and 
Cp |5j. To determine Cf®, Cf®, Cg® and Cf® precisely, we need to measure many observables 
in several radiative and rare semileptonic decays, which potentially can probe new physics at a 
scale of a few TeV. 
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Figure 2. Examples of new physics processes via a charged Higgs boson (left), charginos (middle 
left), neutralinos (middle right) and via Majorana-type neutrino interactions (right). 


Lepton-number-violating decays are highly suppressed in the SM and may need new physics 
processes. Figure [ 2 ] (right) shows a W annihilation diagram into in which the neutrino 
mixes into an antineutrino producing like-sign leptons that are forbidden in SM interactions. 
Such processes require Majorana-type neutrinos that are absent in the SM [ 6 |. 

The decays B° coco and B° — > ucj) also involve FCNC processes that are mediated by 
gluonic penguin loops included in H e g (see Eqn. ([!])). Here, the short-distance contributions are 
parameterized by the Wilson coefficients C 3 , C 4 , C 5 and C%, while the long-distance contributions 
involve the operators O 3 , O 4 , O 5 and 0$. Figure [3] shows the lowest-order diagrams for these 
decays. New physics loops depicted in Fig. [2] may also contribute here. These charmless vector 
vector decays involve three amplitudes. In the transversity frame, these are the longitudinal 
amplitude Aq (S-wave), the transverse amplitude At (P-wave) and the parallel amplitude A\\ 
(D-wave). For measuring CP violation, they need to be known. 



i UJ d d UJ 


Figure 3. Lowest-order diagrams for the B° —> ww color-suppressed tree (left), singlet penguin 
(middle left) and gluonic penguin (middle right) and the gluonic penguin for B° —»• uj(j) (right). 

In chapter [2j we present new BABAR measurements of the direct CP asymmetry in B —> X s 7 
using a semi-inclusive analysis. We extract the ratio of Wilson coefficients Xm(C'| ff /C| ff ) from 
a measurement of the difference in CP asymmetries between charged and neutral B decays. 
We also show CP asymmetry measurements for B — > X s ^7 decays. In chapter [3j we present 
our branching fraction and CP asymmetry measurements of B —> X s £ + £~ decays using a semi- 
inclusive analysis. In chapterMj we summarize our branching fraction upper limits onB-> tt£ + £~ 
and B — > r]£ + £~ . In chapter [5[we summarize our results on searches for lepton number violation 
in exclusive B + — > X~£ + £' + modes. In chapter [ 6 ] we present our results on the charmless vector 
vector decays B° —> ujuj and B° — > ujcj) and in chapter [t] we end with concluding remarks. Note 
that BABAR performs all analyses blinded meaning that results are sealed until selection criteria 
and fitting procedures are finalized. 


2. Measurement of CP Violation in B —> X s ^ 

In the SM, the B — > X s ^ branching fraction is calculated at next-to-next-to-leading order (up 
to four loops) yielding B(B — »■ X s ^) = (3.14 ±0.22) x 10 -4 for photon energies E* > 1.6 GeV in 
the center-oLmass (CM) frame 17;,(8j . For larger values of E* : the prediction depends the shape 
of the of the E* spectrum, which is modeled in terms of a shape function |9| that depends on the 
Fermi motion of the b quark inside the B meson and thus on the b quark mass. Since the shape 
function is expected to be similar to that determining the lepton-energy spectrum in B — >• X u £v, 
precision measurements of the E* spectrum help to determine |I4b| more precisely 
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121. The 


measurement of B(B —> Xg'y) also provides constraints on the charged Higgs mass 0 ! 
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Experimentally, the challenge consists of extracting B —> X s ^y signal photons from those of 
7T° and rj decays, copiously produced in qq continuum (with q = u, d, s , c) and BB processes 
that increase exponentially at smaller photon energies. One strategy consists of summing 
b —> s'} exclusive final states. In a sample of 471 x 10 -6 BB events collected with the BABAR 
at the PEP-II asymmetric storage ring at the SLAC National Laboratory, we 


detector 115,16 


reconstruct 38 exclusive final states containing one or three kaons with at most one Kg, up to 
four pions with at most two 7r°s and up to one ??. We require photon energies in the CM frame 
of 1.6 < E* < 3.0 GeV. Previously, we published total and partial branching fractions 


Here, we focus on the measurement of direct CP asymmetry, which is defined by 


17 


W7) - S(S 


Xsi) - B(B -7 X sl ) 


For this analysis 


18 


A CP (B x s7 ) + B(B 

, we select 16 self-tagging modes, ten B + 


Xsl) 


(2) 


^ + oQ an d six B° —> h °7 hnal 
state^J We maximize the signal extraction using a bagged decision tree with six input variables. 
This improves the efficiency considerably with respect to the standard A E = E* B — E^ eam 
selection where E? and E* B are the beam energy and B meson energy in the CM frame, 
respectively. To remove continuum background, we train a separate bagged decision tree using 
event shape variables. For each X s mass bin, we optimize the sensitivity S/\J(S + B) where 
S(B) is the signal (background) yield using loosely identified pions and kaons. To extract Acp , 

we fit the beam-energy-constrained mass rags = 


E * 2 — 

beam 


p* B simultaneously for H-tagged 


and H-tagged events where p* B is the B momentum in the CM frame. After correcting the raw 
Acp for detector bias determined from the tyies sideband below the signal region, we measure 


A C p(B -7 X sl ) = (1.73 ± 1.93 stat ± lA2 sys )% 
of —0.6% <Acp< 2.8% at 95% confidence level 


181, which agrees well with the SM prediction 
CL) (T9 and which supersedes the old BABAR 
measurement |20|. Though this result is the most precise single direct Acp measurement, the 
uncertainty is sufficiently large to allow for new physics contributions in Cf ff . Figure [ 4 ] (bottom 

The CP asymmetry 
A+7) - A C p(B° -7 


part) shows our result 18 


in comparison to the Belle measurement 1211 
difference between B + and B° decays, A Acp{B —> X s ^) = Acp(B + — 

Xgj), is very sensitive to new physics since it originates from the interference between the 
electromagnetic and the chromomagnetic penguin diagrams in which the latter enters through 
higher-order corrections. Calculations yield [ 19 


a f a (~< e f f 

A Acp(B -7 Xg'y) ~ 47 r 2 a s — -Im 8 ff ~ 0.12- , T Tm s ff , (3) 

v u m b C e 7 ff 100 MeV (* ff 

where A 78 is the hadronic matrix element of the O 7 — Os interference, predicted to he in the 
range 17 MeV < A 78 < 190 MeV. In the SM, A Acp{B -7 V s7 ) vanishes since Cj S and Cf ff 
are real. However in new physics models, these Wilson coefficients may have imaginary parts 
yielding non-vanishing A Acp(B -7- V s7 ) (22-24). From a simultaneous fit to B + and B° modes, 
we measure A Acp(B -7 X s ^)= (5.0 ± 3.9 s t a t ± 1.5 sys )% from which we obtain the constraint 
-1.64 < lm(Cg S /Cj 3 ) < 6.52 at 90% CL. This is the first A Acp(B -7- X s j) measurement and 
first constraint on Zm(Cf ff /Cf ff ). 

Figure 0 (left) shows the Ay 2 of the fit as a function of Zm{C^/C^). The Ay 2 
dependence on Im(Cg /C 7 ) is not parabolic indicating that the likelihood has a non-Gaussian 
shape. The reason is that Ay 2 is determined from all possible values of A 78 . In the region 
0.2 <Zm(C| ff /Cf ff )< 2.6, a change in Zm(C| ff /C' 7 ff ) can be compensated by a change in A 78 
leaving Ay 2 unchanged. For positive values larger (smaller) than 2.6 (0.2), Ay 2 increases slowly 


1 

2 


h + = Kgiv + , A' + 7t°, A' + 7r + 7r , Kgn + n°, K + n°n °, Kgn + n n + , K + iv + n n°, 
h° = K + n~, K + n~n°, K + n + n~Tv~, K + n~Tr°n°, K + n~r l /, K + K + K~tt~ . 


K°gTr + n°n 0 , K+r), K+ K+K~. 
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SM A(B-X std y) 



Figure 4. Summary of Acp measurements for B 
part) from BABAR 18 and Belle 

Belle 


part) from BABAR 
B 


25 


27 


Xs'y 19 and for B —> X s/ pf 


21 

28! 


22 


X^ from semi-inclusive analyses (bottom 
and for B —> X s +d 1 from fully inclusive analyses (top 
and CLEO )29| in comparison to the SM prediction for 
respectively. 


30 



Figure 5. The Ax 2 (left) and A 78 (right) dependence onIm(C§^/C^). The blue dark-shaded 
(orange light-shaded) region shows the 68 % (90%) CL interval. 


(rapidly) since A 78 remains nearly constant at the minimum value (increases rapidly). For 
negative Xm(C| ff /C| ff ) values, A 78 starts to decrease again, which leads to a change in the A % 2 
shape. Figure [5j (right) shows A 78 as a function of Xm(Cf® /C 7 ). 

In the fully inclusive analysis, Acp involves contributions from B —> and B —> X^ that 

cannot be separated on an event-by-event basis. Therefore, we define Acp here as 























































































































(4) 


Acp{B —> X s+d 7 ) 


B{B -+ X s+dl ) - B(B ->• X fl+d7 ) 
B(B -a A s+d7 ) + B(£ ->• X s+d7 )' 


We tag the flavor of the non-signal B flavor by the lepton charge in semileptonic decays 
a sample 384 x 10 6 BB events, we measure Acp{B —> -Xs+d 7 ) = 0.057 ± 0.06^^* ± 0.018 

Figure [4] (top part) shows all Acp{B 


Using 

after 


correcting for charge bias and mistagging 
measurements from BABAR 
prediction (22,30 


27 , Belle 


17 


28 


and CLEO 29 


sys 

Xs+dl ) 

which all agree well with the SM 


3. Study of B — y X s £ + £~ Decays 

Using a senri-inclusive approach, we have updated the partial and total branching fraction 
measurements of B —> X s £ + K modes (I = e or fi) with the full BABAR data sample of 
471 x 10 6 BB events. We reconstruct 20 exclusive final states: K + , Kg, K* + (K + tt° , Kgir + ), 


K*°(K + tt~, KgTT 0 ), K + 7T~ir°, Kgir + Tr~ with Kg -+ 7 r + 7 r” recoiling against e + e~ or 1311. 

After accounting for K\ modes, Kg —> and n 0 Dalitz decays, the selected decay modes 

represent 70% of the inclusive rate for hadronic masses mx s < 1.8 GeV/c 2 . Using JETSET 
fragmentation 


(321 and theory predictions (33-38 , we extrapolate for the missing modes and 
those with mx s > 1.8 GeV/c 2 . We impose the requirements tues > 5.225 GeV/c 2 and 
0.1 (0.05) < A E < 0.05 (0.05) GeV for X s e + e~ (X s r + r,~) modes. We define six bins of 
the momentum-squared transferred to the dilepton system q 2 = m 2 u and four bins in hadronic 
mass mx s • Table [l] shows the defined ranges of these bins. 


Table 1. Definition of the q 2 , mu and mx s bins. 


q 2 bin 

q 2 range [GeV 2 /c 4 ] 

mu range [GeV/c 2 ] 

mx s bin 

mx s range [GeV/c 2 ] 

0 

1.0 - 6.0 

1.00 - 2.45 



1 

0 . 1 - 2.0 

0.32- 1.41 

1 

0.4 - 0.6 

2 

2.0-4.3 

1.41-2.07 

2 

0.6 - 1.0 

3 

4.3-8.1 

2.07 -2.6 

3 

1.0 -1.4 

4 

10.1 -12.9 

3.18-3.59 

4 

1.4 - 1.8 

5 

14.2 - {jur — m * K ) 2 

3.77 - {tub — pi*k) 




To suppress e + e _ —> qq and BB combinatorial background, we define boosted decision trees 
(BDT) for each q 2 bin separately for e + e - and /i + fj,~ modes. From these BDTs, we determine a 
likelihood ratio (Lr) to separate signal from qq and BB backgrounds. We veto J/i/> and ^{2S) 
mass regions and use them as control samples. We measure dJ3(B -+ X s £ + l~/dq 2 in six bins of 
q 2 and four bins of mx 3 ■ We extract the signal in each bin from a two-dimensional fit to mES 
and Lr. As examples, Figs, [h] and [7] show the mES and Lr distributions for e + e _ modes in bin 
g 5 and for fi + ^ modes in bin q\, respectively. Clear signals are visible both in the rngg and 
Lr distributions. Figure | 8 | shows the differential branching faction as a function of q 2 (left) and 
mx s (right) [31 . Table [2] summarizes the differential branching fractions in the low and high 
q 2 regions in comparison to the SM predictions (34,38-47]. In both q 2 regions, the differential 
branching fractions are in good agreement with the SM predictions. These results supersede the 
previous BABAR measurements 


48 and agree well with the measurements from Belle 49 


The direct CP asymmetry is defined by 
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Figure 6. Distributions of tues (left) and likelihood ratio (right) for B — > X s e + e~ in q 2 
bin gs showing data (points with error bars), the total fit (thick solid blue curves), signal 
component (red peaking curves), signal cross feed (cyan curves), BB background (magenta 
curve), e + e _ —> qq background (green curves) and charmonium background (yellow curves). 




Figure 7. Distributions of rriES (left) and likelihood ratio (right) for B — > X s qi + [i~ in q 2 
bin q\ showing data (points with error bars), the total fit (thick solid blue curves), signal 
component (red peaking curves), signal cross feed (cyan curves), BB background (magenta 
curve), e + e _ —> qq background (green curves) and charmonium background (yellow curves). 



Figure 8. Differential branching fraction of B X s e + e~ (blue points), B — > X s [i + /r~ (black 
squares) and B —> X s i + t~ (red triangles) versus q 2 (top) and versus mx s (bottom) in comparison 
to the SM prediction (histogram). Grey-shaded bands show J/-0 and ^{2S) vetoed regions. 
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Table 2. Measured B —> X s £ + £~ branching fractions in the low and high q 2 regions from 
BABAR 1 31 1 and the SM predictions. Uncertainties are statistical, systematic and from model 
dependence, respectively. 


Mode 

BABAR [10 -6 ] 

SM [10- 6 ] 

BABAR [10 -6 ] 

SM [10- 6 ] 

q 2 [ GeV 2 /c 4 ] 

1 - 6 

1 - 6 

> 14.2 

> 14.2 

B -> X s p + pL~ 
B -7 X s e + e~ 

B -7 X s £ + £~ 

0.66±g;f|±g;i2 ± 0.07 

1 - 93 ^ 45 ™ ±0.18 
l-60i^; 1 1 3 7 ±0.07 

1.59 ±0.11 

1.64 ±0.11 

0 fin +°- 31 + 0 - 05 4- 0 (i(i 
U.UU 0 . 29 — 0.04 ^ u ' uu 

0 cm+°- 19 + 0 - 03 4- 0 0(1 
U . 00 ^ 0 . 18 - 0.03 ^ U - UU 

0 ^7+°- 16 +°- 03 4-0 0(1 
u - 0/ - 0 . 15 - 0.02 111 u ' uu 

0 25 -1-0 ' 07 
u - zo - 0.06 

0 25 -1-0 ' 07 
u - zo - 0.06 

0 25" 1-0,07 
u - zo - 0.06 




0.5- 


-0.5 — 


J/\|/ 


10 


15 


20 

q 2 (GeV/c 2 ) 


25 

, 2\2 


Figure 9. The CP asymmetry versus q 2 . Grey bands show the J/i[) and i/j(2S) vetoed regions. 


We use 14 self-tagging modes consisting of all B + modes and B° modes with decays to a 
K + to measure Acp{B — > X s £ + £~) in five q 2 bins. Due to low statistics, we have combined 
bins (/4 and q§. Figure [ 9 ] shows the CP asymmetry as a function of q 2 . The SM prediction of 


54 


55 


50 


53 


the CP asymmetry in the entire q 2 region is close to zero 
however, Acp may be significantly enhanced 
A CP (B -7 X s £+£~) = 0.04 ± 0.11 stat ± 0.01, f _ 

prediction. The CP asymmetries in the five q 2 bins are also consistent with zero. 


In new physics models, 
In the full range of q 2 , 


we measure 


31 , which is in good agreement with the SM 


4. Search for B —» -n/■q£ + £~ Decays 

In the SM in lowest order, B —> Xd£ + £~ modes are also mediated by the electromagnetic penguin, 
Z penguin and WW box diagrams. However, they are suppressed by \Vtd/Vt s \ 2 — 0.04 with 
respect to the corresponding B -7- X s £ + £~ decays. In extensions of the SM, rates may increase 
significantly 56 . Using 471 x 10 6 BB events, we recently updated the search for B —> n£ + £~ 
modes and performed the first search for B —> rj£ + £~ modes [57j. The SM predictions he in the 
range B(B -> n£ + £~) = (1.96-3.30) x 10 -8 and B(B —> rj£ + £~) = (2.5-3.7) x 10 -8 where the 
large ranges result from uncertainties in the B —>• it form factor calculations [56,58,59] and from 
a lack of knowledge of B —> q form factors 
We fully reconstruct four B 
77 and 7] —> 77 , 7 r + 7 r _ 7 r° recoiling against e + e _ or /r + // _ . We select leptons with pi > 0.3 GeV/c, 
recover eX bremsstrahlung losses, remove 7 conversions and require good particle identification 


respectively. 

■n£ + £~ and four B —> rj£ + £~ final states by selecting tt^ , 7r° 


















































for eF , and vr^. We select photons with E~ f > 50 MeV and impose 7 r° and q mass constraints 
of 115 < m 77 < 150 MeV/c 2 and 500 (535) < m 77 {m^) < 575 (565) MeV/c 2 , respectively. For 
the ?7 —>• 77 final state, we require (£i >7 — -E 2 i 7 )/(F % 7 + F? 2 , 7 ) < 0.8 to remove asymmetric 
qq background that peaks near one. We veto J/%b and tp(2S) mass regions and use four 
neural networks (NN) to suppress combinatorial BB and qq continuum backgrounds, separately 
for e + e~ and for /x + /i _ modes. The NNs for suppressing BB background use 15 (14) input 
distributions for e + e _ (// + /x“) modes, while those for suppressing qq continuum use 16 input 
distributions for both modes. For validations of the fitting procedure and peaking backgrounds, 
we use pseudo-experiments and the vetoed J/tp and ip(2S) samples. 

For B + —> 7 t + £ + £~ and B° —> tt°£ + £~ , we perform simultaneous unbinned maximum 
likelihood (ML) fits to tues and A E distributions for e + e - and /u + /i - modes separately. We 
include the B + —> K + £ + £~ mode in the fit to extract the peaking background contribution in 
the B + —> tt + £ + £~ modes by reconstructing the K + as a 7 r + . For B° —> r ]£ + £~, we perform 
simultaneous unbinned ML fits to tues and A E distributions, again for e + e _ and n + modes 
separately. In addition, we perform fits for the isospin-averaged modes B —> 7 re + e _ and 
B —> 7 lepton-flavor-averaged modes B + —> tt + £ + £~, B° —> tt°£ + £~ and B° —> r]£ + £~ 
and both isospin- and lepton-flavor-averaged modes B —> tt£ + £~ . 



Branching fraction [1 O' 8 ] 


Figure 10. Branching fraction upper limits at 90% CL for B —>• n£ + £ 
from BABAR 157 and Belle 61 and the measurement of B + —> ir + £ + £~ 


and B —> ij£ + £ 
from LHCb 


modes 


We see no signals in any of these modes and set branching fraction upper limits at 90% 
CL. Figure [Tol shows them in comparison to results from Belle 1611 and a measurement of 
B{B+ -> tt+TV) from LHCb 6§. For B° -> n°£ + £- our branching fraction upper limit is 
the lowest and so far only BABAR has searched for B° —> ij£ + £~ modes. The present branching 
fraction upper limits lie within a factor of two to three of the SM predictions. 


5. Search for Lepton Number Violation in B + —> X~£ + £' + Decays 

In the SM, lepton number is conserved in low-energy collisions. However, in high-energy and 
high-density interactions, lepton number may be violated |63). Many models beyond the SM 
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predict lepton number violation (LNV) with rates j64| that may be accessible already in present 
data samples. These models also predict Majorana-type neutrinos |6j for which particles and 
antiparticles are identical. Via oscillation of a neutrino into an antineutrino, lepton-number 
violating decays become possible such as B + —> K~£ + £'^ depicted in Fig. [ 2 ] (right). The 
observation of atmospheric neutrino oscillations confirms hat neutrinos carry mass 65] but we 
do not know if any Majorana-type neutrinos exist. However, lepton number violation is also a 
necessary condition to explain the observed baryon asymmetry in the universe [66]. 

Using the full BABAR data set of 471 x 10 6 BB events collected at the T(45') peak, we 
have searched for lepton number violation in 11 B + decay^] |67 . We select events with 
more than three charged tracks of which two are identified as like-sign leptons having a 
combined momentum less than 2.5 GeV/c in the laboratory frame. We remove e + and e~ 
from photon conversions. We define sufficiently wide mass regions around the p~ (0.47 < 
m V7r o < 1.07 GeV/c 2 ), K*~ (0.770 < < 1.01 GeV/c 2 ) and D~ (1.835 < mKir-K < 

1.895 GeV/c 2 ) mesons to allow reasonable modeling of backgrounds. We combine the 
tt ~, K~ , 7 r _ 7 r°, K~n°, and ~ candidates with the two leptons to form a B 

candidate. We remove l±hA combinations with an invariant mass close to that of the J as 
the hA may be a misidentihed pA from — >• p + p~ . The misidentification rate is about 2 %. 


Table 3. Summary of signal yield after fit bias correction with statistical uncertainty, 
reconstruction efficiency e, significance S including systematic error, measured branching fraction 
B x 10' and the 90% CL upper limit Bjjl x 10 7 . Yields and efficiencies for the B + -» K*~£ + £' + 
modes are for K~ir° and Y^ 7 r _ final states, respectively. 


Mode 

Yield 

e 

S 

B 

Bui 


[events] 

[%] 

M 

[ 10 - 7 ] 

[ 10 - 7 ] 

K*~e + e + 

3.8 ±3.3, 0.8 ±3.9 

11.5 ±0.1, 12.1 ±0.1 

1.2 

1.7 ± 1.4 ±0.1 

4.0 

K*~e + p + 

-1.9 ±4.7, -5.1 ±2.6 

7.9 ±0.1, 8.5 ±0.1 

0.0 

-4.5 ±2.6 ±0.4 

3.0 

K*~p + p + 

2.3 ±1.8, 2.0 ±1.8 

6.1 ±0.1, 5.8 ±0.1 

1.3 

2.4 ±1.8 ±0.4 

5.9 

p~ e + e + 

-2.1 ±5.7 

12.1 ± 0.1 

0.0 

-0.4 ± 1.0 ±0.1 

1.7 

p~e + p + 

4.6 ± 11.4 

10.3 ±0.1 

0.4 

1.0 ±2.4 ±0.2 

4.7 

p~p + p + 

2.9 ± 6.8 

7.3 ±0.1 

0.5 

0.9 ±2.0 ±0.3 

4.2 

D~e + e + 

3.9 ±4.8 

10.2 ± 0.1 

1.0 

8.8 ± 8.6 ± 1.5 

26 

D~e + p + 

1.1 ±3.2 

7.7 ±0.1 

0.5 

3.4 ± 9.4 ± 1.1 

21 

D~p + p + 

—1.7 ± 2.5 

5.7 ±0.1 

0.0 

-6.5 ±9.9 ±0.9 

17 

K~e + p + 

5.5 ±3.5 

15.2 ±0.1 

1.8 

0.6 ±0.5 ±0.1 

1.6 

7r~e + p + 

3.8 ±3.5 

16.4 ±0.2 

1.2 

0.5 ±0.5 ±0.1 

1.5 


For each signal mode, we construct BDTs to discriminate signal from BB and qq backgrounds 
using nine inputs consisting of event shape variables, kinematic observables, flavor tagging and 
the proper decay time. If more than one candidate is found, we choose the one with the 
smallest x 2 in the fit to the B decay vertex. We perform a simultaneous unbinned ML fit 

3 B+ p~ (tt~ n 0 )£ + £' + , B+ -> K*-(K° s tt- and K~n°), )£+£'+, B+ D~ (K~tv + tt-)£+£'+ , B + n~e+p + 
and B + — > n~e + p + where l + ,£' + = e + or p + . 
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to rriES-, AE and the BDT output distributions. For the B + —> p £ + £' + , B + —> K* £ + £' + 
and B + —> D~£ + £ ,+ final states, we include the m^ir, rriKn and iJiKinr mass distributions, 
respectively. The background PDFs consist of an Argus function [68] for mESi first- or second- 
order polynomials (for K ~, tt~, p ~, K*~ modes) or a CruijfQ function (for D~ modes) for 
A E, a non-parametric kernel estimation KEYS algorithm 1691 for the BDT output and a first- 
order polynomial plus a Gaussian function for the resonance masses. The corresponding signal 
PDFs consist of a Crystal Ball function [70] for tres , a Crystal Ball function plus a first-order 
polynomial (for modes with 7T°s) for A E, the simulated distribution in form of a histogram 
for the BDT output and for the D , K*, p masses two Gaussians, a relativistic Breit-Wigner 
function and a Gounaris-Sakurai function |71 , respectively. We checked the fit procedure with 
a simulated background sample having the same size as the on-resonance data sample. We 
further performed a blinded fit to the on-resonance data sample confirming that the background 
distributions agreed with the background PDFs. Selection efficiencies vary between 6% and 16% 
depending on the final state. 



5.24 5.25 5.26 5.27 5.28 5.29 

m ES (GeV/c?) 




Figure 11. Projections of the multidimensional fit onto hies (left)) AE (middle) and the BDT 
output (right) for B + —> ir~e + p + showing data (points with error bars), total fit (solid blue 
line), signal PDF (green histogram) and background (dashed magenta line). 


Figures 11 and 12 show projections of the fit on the discriminating variables for B + —> 
and B + —> K*~(KgTT~)p + p + , respectively. Table 0 summarizes our results. In all 11 


7r e 

modes, the data are consistent with combinatorial background. We see the highest significance 
of 1.8cr in B + —> K~e + p + . We set Bayesian upper limits on the branching fraction at 90% CL 
using a fiat prior (see Tab. [3]). The additive systematic uncertainty that includes contributions 
from the PDF parameterization, fit biases, background yields and efficiencies is mode dependent 
between 0.2 and 0.7 events. The total multiplicative uncertainty on the branching fraction is 
5% or less. The branching fraction upper limits at 90% CL lie in the range 1.5 - 26 xl0~' 
where the lowest limit is set in the B + —> n~e + p + mode. Figure 13 summarizes all results 
of lepton-number-violating B decays from BABAR 
including results for B° 


A+T, B~ -a A£-. B- 


67] , Belle 

> M 


72 , LHCb 


75 and B 


IT 


73 
/K 


and CLEO [74] 
All 


-£+£+ 76 


limits are set at 90% CL except for LHCb whose limits are set at 95% CL. 


6 . Study of B° -A UJUJ and B° —> ojcj) Decays 

The longitudinal polarization fraction El in charmless vector vector decays poses a puzzle. In 
tree-dominated decays like B° -A p + p~ and B + -A p + p° , El is nearly 100% while in decays 
with dominant penguin contributions like B -A K*p , El is around 50% [77] . Are there large 
transverse SM contributions that reduce El [78|-83] or is this caused by new physics 82,84-94 ? 


Thus, it is interesting to investigate other charmless vector vector decays such as the so far not- 
observed modes B° -A wu and B° -A uucf) 95 . In the SM, the branching fractions are expected 


4 The Cruijff function is a centered Gaussian with different left-right resolutions and non-Gaussian tails: 
f(x) = exp ((x - m) 2 /(2a 2 L R + a Ll n(x - m) 2 ). 
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Figure 12. Projections of the multidimensional fit onto tues (top left), A E (top right), the 
BDT output (bottom left) and K*~ mass (bottom right) for B + —> K*~ (Kg7r~)e + n + showing 
data (points with error bars), total fit (solid blue line), signal PDF (green histogram) and 
background (dashed magenta line). 



Figure 13. Branching fraction upper limits at 90% CL for LNV decays from BABAR [67] (solid 
blue points), Belle [72j (solid red squares) and CLEO [74] (black diamonds). In addition, LHCb 
upper limits at 95% CL [73| (open triangles) are shown. 


to be of the order of 0(10 6 ) for B° —> ww and 0(10 7 ) for B° 
longitudinal polarization fractions of J~l> 0.8 for both modes [93,94 


■> ucj). The SM predicts 
Charmless vector vector 
modes are also well suited to measure the Unitarity Triangle angle a [96,97 . The Scan Method 
group has determined a — j3 contours from a x 2 fit to measured branching fractions, longitudinal 
polarizations and CP asymmetries using all observed charmless vector vector decays 98 . The 


decay amplitudes of each mode are expressed in terms of tree, color-suppressed tree, gluonic 
penguin, singlet penguin, electroweak penguin and kP-annihilation/IT- exchange amplitudes. 
For decays involving K* s, SU(3) breaking is taken into account. All contributions up to order 
0(A 5 ) are considered where A = sin0 c (Cabibbo angle), since the leading amplitude is already 
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at order 0(A 3 ). Figure 14 shows the 90% CL a-/3 contour determined from the fit. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


p 

Figure 14. The a-f3 contour obtained from a ML fit to branching fractions, longitudinal 
polarizations and CP asymmetries measured in all observed B —> VV modes. The magenta 
band shows the result of the sin 2/3 measurements in B —> ccK^ L decays 


77 


we reconstruct the F-daughter 
K+K~ 


99 


If multiple 


Using the full BABAR data sample of 471 x 10 6 BB events, 
candidates via their decays oj —> with 7r° —> 77 and 0 

candidates exist, we select the one for which the \' 2 probability of a fit to the two vector meson 
masses is smallest. The combinatorial background from e + e _ —> qq collisions dominates. We 
use a tight selection on the angle 0t between the thrust of the B signal candidate in the T(4S') 
rest frame and that of the rest of the event, requiring |cos0-r| < 0.8 (0.9) for B° — > loco and 
B° —> u}(j) decay modes. Furthermore, we define a Fisher discriminant T based on four shape 
and kinetic variables as inputs. We perform an extended unbinned ML fit to extract the signal 
and background yields from the data. We define the PDF as a product of six individual PDFs 
including m^s, A E, J 7 , masses and helicity angles of the two vector mesons and the decay angle 
0 between the 7r° in the dipion rest frame and the u flight direction: 

V) = Vjimiss) x Vj(AE l ) x V^F) x X V0^) (6) 

where the last term is not present in B° —> w0. For signal, we use a sum of two Gaussians for 
mES and A E, a two-piece normal distribution for E, relativistic Breit-Wigner functions for the 
my distributions, each convolved with two Gaussians to account for detector resolution. We 
parameterize the helicity angles by the angular distribution: 


W(0y 1 ,0y 2 ) oc (1/4)(1 — J-jj) sin^ Qy sin^ dy 2 + Pj J cos 2 9y 1 cos 2 0y 2 


(7) 


convolved with a resolution function for each angle. Thus for signal, the PDF is factorized as 
Vj {rriy ,, rriy 2 , 9\ , 6\) = Vj i'l'i'iy ) x Vj(my 2 ) X W(9y 1 ,9y 2 ). For the angles ' 01 , 2 ) the PDFs are 
sin 2 2 distributions. For combinatorial background, we use an Argus function 68 for rriESi 


a second-order polynomial for A E and a two-piece normal distribution for T. The masses and 
helicity angles of the two vector mesons are considered to be independent and thus are factorized. 
We use a third-order polynomial for the my distributions. The PDFs of the helicity angles are 
third-order polynomials for combinatorial background where the parameterization is obtained 
from on-peak sideband data (tues < 5.27 GeV/c 2 ). Similarly, we use a third-order polynomial 
in cos 0 for the u> decay angle. We also consider BB peaking background determined from 
simulation. We parameterize mESi AF and T with similar functions as those for the signal. The 
PDFs for my and 0 angle are similar to those of the combinatorial background. For the helicity 
angles, we use a fourth-order polynomial. For the signal, BB and qq background components, 
we determine the PDF parameters from simulation. We study large control samples of B —> Dn 
decays of similar topology to verify the simulated resolutions in A E and vies , adjusting the 
PDFs to account for any differences found. 
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Figures [15] (left, middle) show the tues and A FI projections of the multidimensional ML fit for 
B° —> lulu (top) and B° — > ljc/) (bottom. Figure |l5| (right) shows the corresponding distributions 


of —2 ln(C(B) / Co) for zero signal normalized to the likelihood at the minimum value for the two 
modes. We fix T to 0.88 but vary it from 0.58-1.0 in the systematic error determination. Signal 
and background yields are extracted from the fit. For B° —> lulu, we observe = 69^J®'2 
signal events. Including systematic uncertainties, this yields a 4.4cr effect. For B° —> Luq i, we 
see no signal. The fit yields = —2.8^4 q signal events. Besides the uncertainty from the 
J~L variation, systematic uncertainties include contributions from the PDF parameterization, 
selection efficiency, use of control samples, number of B mesons and B-daughter branching 
fractions. The largest systematic error comes from yield bias estimation in the fit yielding order 
0(5 events). We measure a branching fraction of B(B — > lulu) = (1.2 ± 0.3^q'|) x 10 -6 , which 
agrees well with the SM prediction. We do not have enough data to measure J~l in this mode. 
For the B° —> luc/) mode, we set a Bayesian upper limit of B(B —> ucj)) = 0.7 x 10 -6 at 90% CL 
assuming a flat prior. 



m es 





rriEs 




Figure 15. Fit results for B° —> lulu (top row) and B° — > cu</> (bottom row) showing rriES (left) 
and A E (middle) projections of the ML fit for data (points with error bars), the total PDF 
(solid blue line), background PDF (black dashed line) and signal contribution (red dashed line). 
The distributions of — 2 ln(C(B)/Co) versus branching fraction (right) for statistical errors only 
(dashed red curve) and statistical and systematic errors combined (solid blue curve). 


7. Conclusions 

We have performed the most precise direct CP asymmetry measurement in the semi-inclusive 
B —»• A’ s 7 decay. The measurement error is still suffiiciently large allowing to accommodate 
new physics contributions in the Wilson coefficient Cf ff . We performed the first determination 
of Xm(C'| ff /C'f ff ) by measuring the difference in CP asymmetries between charged and neutral 
B decays. We measured partial branching fractions and CP asymmetries in B —> X s £ + £~ that 
agree well with the SM predictions. We find no evidence for B —> tt/'ij£ + £~ decays and set 
branching fraction upper limits at 90% CL that are a factor of two to three above the SM 
predictions. We also find no evidence for lepton number violation in B + — > X~£ + C + decays and 
set stringent branching fraction upper limits at 90% CL. We find first evidence for B° — > lulu 
decays and measure a branching fraction that is consistent with the SM prediction. However, the 









14 


data sample is too small to extract J~l from an angular analysis. We set an improved branching 
fraction upper limit at 90% CL for B° —> ojcj). Significant improvement on these measurements 
are expected to come from Belle II and for some decays also from LHCb. 
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